Giant neutron halos in the non-relativistic mean field approach 
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Giant neutron halos in medium-heavy nuclei are studied in the framework of the Hartree-Fock- 
Bogoliubov (HFB) approach with Skyrme interactions. The appearance of such structures depends 
sensitively on the effective interaction adopted. This is illustrated by comparing the predictions of 
SLy4 and Skl4 in the Ca and Zr isotopic chains. The former force gives no halo effect, the latter 
predicts a neutron halo in the Zr chain with A>122 due to the weakly bound orbitals 3pl/2 and 
3p3/2. The structure of the halo is analyzed in terms of the occupation probabilities of these orbitals 
and their partial contributions to the neutron density. The anti-halo effect in Ni and Zr isotopes 
is also discussed by comparing the occupation probabilities of Hartree-Fock neutron single-particle 
states near the Fermi energy with the corresponding HFB quasiparticle states. 

PACS numbers: 21.10.Gv, 21.10.Pc, 21.60.Jz, 27.40,+z, 27.50.+C, 27.60.+j 
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I. INTRODUCTION 

In the present years several projects for the construc- 
tion of a new generation of radioactive beam facilities are 
in progress (see, e.g.,Ref. 0|). Such facilities will permit 
to investigate the properties of unstable nuclei situated 
close to the drip line regions. From the theoretical side, 
many efforts are devoted to perform accurate predictions 
to locate the proton and neutron drip lines as well as to 
describe the behavior of unstable nuclei. Unfortunately, 
theoretical predictions in exotic regions of the nuclear 
chart can be rather model dependent. 

Self-consistent mean field methods are well suited 
theoretical tools for describing medium and heavy nu- 
clei. There are two main lines of investigations based 
on the mean field approach, namely the relativistic 
mean field (RMF) method where effective Lagrangians 
arc treated in the tree approximation, and the non- 
relativistic Hartree-Fock (HF) method using effective in- 
teractions like the Gogny and Skyrme forces. Recent re- 
views can be found in Refs.0, 13 ■ When approaching 
the drip lines, one deals with open shell nuclei where 
the effects of pairing correlations become quite impor- 
tant, especially for such properties as the tails of matter 
distributions. The pairing correlations can be described 
by the non-relativistic Hartree-Fock-Bogoliubov (HFB) 
theory 0, |f|, or by the relativistic Hartree-Bogoliubov 
(RHB) theory @, and one must add to the models a phc- 
nomenological pairing interaction acting in the particle- 
particle channel. 

Furthermore, the chemical potential A becomes close to 
zero in the vicinity of drip lines and it is necessary to treat 
properly the contributions of the quasiparticle continuum 
when evaluating the pairing correlations 0. Thus, the 
most appropriate approach for such cases is to solve the 
self-consistent mean field equations in coordinate space, 
and this is the method we use to obtain the results of 



this work. 

Important discrepancies are often found in the position 
of the neutron drip line predicted by different models. 
For instance, these differences clearly appear in the neu- 
tron drip line of Ni isotopes calculated in relativistic 
and non-relativistic^ approaches. Even among different 
paramctrizations of effective Skyrme interactions used in 
HFB approach one finds different drip line predictions, as 
illustrated in the recent work of Ref. Q and in the results 
of the present work. While waiting for new experimental 
data which will help to discriminate among the available 
models it can be interesting to analyze the reasons for 
some of these discrepancies. This is one of the aims of 
this paper. 

A very interesting phenomenon has been recently pre- 
dicted within the RHB approach, the formation of a neu- 
tron giant halo (with up to six neutrons involved) in 
some very neutron-rich isotopes. By plotting the neu- 
tron distribution radius with respect to A a kink ap- 
pears where the halo structure starts being formed. This 
effect has been found to be particularly strong in Ca 
(with A > 60) @ and Zr isotopes (with A > 122) [Toj . 
These predictions are based on the NLSH parametriza- 
tion, and similar results are obtained with the TM1 
parametrization^lj- The giant halo phenomenon is also 
found in the near-drip line Zr isotopes if one uses the NL3 
parametrization in an RMF plus resonant continuum 
BCS0I- Until recently, the giant halo effect in medium- 
heavy nuclei has been much less investigated within the 
non-relativistic mean field approach, apart from Ref.Q 
where the study of halo is concentrated on the Ca iso- 
topes. In this paper, we focus our discussion on the Zr 
isotopes where we find that, in the HFB approach with 
Skyrme forces the position of the neutron drip line de- 
pends on the particular force adopted and consequently, 
the giant halo effect seems to be more model dependent 
than in relativistic approaches. 
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A well known phenomenon for which relativistic and 
non-relativistic predictions are quite different is the 
charge isotope shift in Pb isotopes: the kink in the charge 
radius at the shell closure N = 126 is not reproduced by 
non-relativistic calculations with the Gogny interaction 
|l3j | or with standard Skyrme parametrizations while it 
can be well reproduced by RMF calculations. Reinhard 
and Flocard |l4j have analyzed this discrepancy with 
a special attention to the spin-orbit paramctrization of 
Skyrme interactions. New parametrizations have been 
proposed in order to reproduce the observed kink of the 
charge isotope shift. If one introduces two parameters 
in the spin-orbit part of the energy functional instead of 
one as in the standard Skyrme interactions, it is possible 
to obtain a spin-orbit potential with isovector proper- 
ties similar to that of RMF. Then, a fitting procedure of 
all the Skyrme force parameters is performed to produce 
new parametrizations which can describe well a chosen 
set of nuclei and also reproduce the kink in the Pb iso- 
tope shift. In the present calculations we have selected 
the paramctrization SkI4, and for comparison purposes 
we also use the more standard set SLy4. 

In this work we investigate the giant halo effect in 
Ca and Zr isotopes within the non-relativistic Skyrme- 
Hartrce-Fock-Bogoliubov mean field approach. We find 
that the neutron drip line of Zr isotopes calculated with 
SkI4 is far enough to give rise to a giant halo phe- 
nomenon. On the other hand, neither SLy4 nor SkI4 can 
produce bound Ca isotopes which are enough neutron- 
rich to lead to a halo effect. The anti-halo effect due to 
pairing correlations is discussed for Zr and Ni iso- 
topes by comparing HFB results with the corresponding 
HF ones. 

The article is organized as follows. In Sec. II we 
briefly describe the theoretical framework. In Sec. Ill we 
present the calculated two-neutron separation energies, 
neutron radii and single-particle spectra for neutron-rich 
Ca and Zr isotopes and we compare our non-relativistic 
results with the correspond ing RH B results obtained with 
the NLSH paramctrization [glllCj. In Sec. Ill we consider 
in more detail the case of Zr isotopes and investigate giant 
halo and anti-halo effects by analyzing occupation proba- 
bilities, neutron density profiles and radii. The anti-halo 
effect is also studied for Ni isotopes. In Sec. IV our 
conclusions are drawn. 



II. THEORETICAL FRAMEWORK 

The theoretical framework used in this paper is the 
Hartree-Fock-Bogoliubov (HFB) approach. For local 
two-body forces and spherical symmetry the HFB equa- 
tions have the following form: 



[h(r) - X]u(E, r) + A(r)v(E, r) = Eu{E, r) 
A(r)u(E, r) - [h(r) - X)v(E, r) = Ev(E, r), 



(1) 



pairing potential; u and v are the upper and lower compo- 
nents of the quasiparticle wave function associated with 
the energy E. 

The quasiparticle spectrum contains a discrete part 
for E less than —A and a continuous part for E above 
—A. To calculate the continuum spectrum the HFB 
equations should be solved with scattering type bound- 
ary conditions for the upper components of the HFB 
wave functions ||. Since the continuum- HFB calcula- 
tions are rather heavy, usually the continuum spectrum 
is discretized by imposing box boundary conditions, i.e., 
the condition that the HFB wave functions vanish at a 
given distance from the nucleus. This is sufficient for our 
present purpose provided that the box radius is properly 
chosen. We have checked that our results obtained with 
a box radius of 20 fm are very close to those of a full 
continuum calculation. 

In the HFB calculations presented in this paper the 
mean field is calculated with a Skyrme type force while 
for the pairing channel we use a zero range interaction 
with the following density dependence: 



V(T! - r 2 ) = V a 



1-x 



P( r ) 
Po 



S(ri-r 2 ). (2) 



where A is the Fermi energy, h is the sum of the kinetic 
energy and the HF mean field potential, and A is the 



We have adopted the values po=0.16 fm -3 , a;=0.5, 7=1. 
The strength Vb is chosen so as to reproduce the gaps 
extracted from the odd-even mass differences ( in the 
regions where such experimental data are available). 

As for the Skyrme interaction, we use two forces, 
i.e., SLy4 and SkI4. The force SLy4, which includes 
constraints coming from the neutron matter equation 
of state, is commonly employed for the description of 
neutron-rich nuclei. The force SkI4 has a spin-orbit po- 
tential similar to the RMF models and it includes in the 
fitting protocol the requirement to reproduce the charge 
isotope shift for Pb isotopes. Due to these characteris- 
tics one expects that this Skyrme force will predict for the 
nuclei close to the drip line properties similar to those of 
relativistic models. We are particularly interested in the 
giant halo which might be formed in heavy nuclei close 
to the neutron drip line. Neutron giant halos were pre- 
dicted by the RHB calculations in Ca and Zr regions for 
A > 60 and A > 122, respectively Q. Thus, the pairing 
strengths Vb are adjusted to be -365 MeV ■ fm 3 and -290 
MeV • fm 3 (-350 MeV • fm 3 and -300 MeV • fm 3 ) for the 
Ca and Zr isotopes calculated with SLy4 (SkI4). 

The Skyrme model predictions will be analyzed in the 
next section by using the HFB approach. The HFB cal- 
culations are performed in spherical symmetry, like in the 
RHB approach. 
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FIG. 1: Two-neutron separation energies for Ca isotopes cal- 
culated with SLy4-HFB (triangles), SkI4-HFB (circles) and 
RHB (diamonds). The RHB results correspond to Ref. 



III. RESULTS OF HFB CALCULATIONS 

A. Separation energies and neutron radii 

The two-neutron separation energy is defined as 

S 2n (N,Z)=E(N,Z)-E(N-2,Z), (3) 

where E(N, Z) is the total energy of the isotope with 
N neutrons and Z protons. The two-neutron separation 
energies for Ca and Zr isotopes are shown in Figs. 1 
and 2. In these figures are displayed only the separation 
energies for the bound nuclei (i.e., those having a neg- 
ative chemical potential in the HFB calculations). The 
results correspond to the box-HFB calculations, which 
give practically the same S2n values as continuum-HFB 
calculations. 

The most important fact we can observe in Fig. 1 is 
the large difference between the drip line location pre- 
dicted by the Skyrme-HFB with SLy4 or SkI4, and the 
RHB calculations of Ref.Q for the Ca isotopes. It can 
be seen that, for both Skyrme forces the drip line is lo- 
cated at 62 Ca, while for RHB calculations the drip line 
extends up to 72 Ca. However, it is found in Ref.jg] that 
the two-neutron separation energy is still positive in 78 Ca 
with the interaction SkM*. This strong sensitivity of the 
drip line location to the particular Skyrme force makes 
difficult the study of neutron-rich Ca isotopes. In the 
Zr isotopes there are also differences among the model 
predictions. As seen in Fig. 2, for the force SLy4 the drip 
line is located at 122 Zr while for SkM at 138 Zr. The lat- 
ter result is similar to the RHB one where the drip line 
is located at 140 Zr. 

Next we analyze the neutron radii, shown in Fig. 3. As 
expected, for Ca isotopes the HFB calculations with the 
interactions SLy4 and SKM do not predict a giant halo 
because the neutron drip line is reached before the halo 



FIG. 2: Two-neutron separation energies for Zr isotopes. The 
symbols are the same as in Fig. 1. The RHFB results corre- 
spond to Ref. 01 . 
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FIG. 3: Neutron radii for Ca and Zr isotopes. The symbols 
are the same as in Fig. 1 and the RHB results correspond to 

Ref. sua. Radii calculated with the formula roA 1 ^ 3 (ro = 1 

fm) are shown by a dotted line. 



structure starts to be formed. The same is happening 
for Zr isotopes if one uses the SLy4 force. However, for 
the force SkM the situation is very different. Thus, as 
seen in Fig. 3, the neutron radii given by the SkM- HFB 
calculations show a change of slope around N = 82 and 
then follow the RHB results. 

To understand better the behaviour of the radii, in 
Figs. 4 and 5 are shown the HF energies of the bound 
states close to the Fermi level. For Ca isotopes we can 
see that, apart from a different splitting of the 1/ states, 
the two Skyrme forces give a rather similar structure for 
the bound spectrum. In both cases the state 1^9/2 is not 
bound, at variance with the RMF calculations in which 
this state becomes weakly bound at A = 62 Due to 
this fact the drip line in the RHB calculations is extended 
up to the region where the giant halo can be formed, in 
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FIG. 4: SkI4-HF (left) and SLy4-HF (right) neutron single 
particle energies for bound states in Ca isotopes. 



FIG. 6: Occupation profiles W(E) (in units of MeV" 1 ) for 
the states 3pl/2, 3p3/2 and 2/7/2 in 132 Zr, calculated with 
SkI4-HFB. 
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FIG. 5: SkI4-HF (left) and SLy4-HF (right) neutron single 
particle energies for bound states in Zr isotopes. 



contrast with the HFB calculations of Ca isotopes based 
on SLy4 or SkI4. It must be noted, however, that this 
tendency is not obeyed by the SkM* parametrization 
which gives a bound l.g9/2 orbital in the Ca isotopes 
with A>60 resulting in a drip line located at higher Ay|. 

For Zr isotopes the structure of the bound spectrum is 
not the same for the two Skyrme forces. We can observe 
that, for the force SLy4 there is only one bound state 
close to the Fermi level, i.e., lftll/2, while for the force 
SkM the states 2/7/2, 3p3/2 and 3pl/2 appear as bound 
states as well. These three weakly bound states have a 
similar structure in the RMF calculations. This explains 
why the SkI4-HFB and the RHB calculations provide a 
similar halo in Zr isotopes. How these weakly bound HF 
states contribute to the halo structure in the presence of 
pairing correlations is discussed in the next subsection. 



B. The structure of the giant halo 

In the presence of pairing correlations the bound 
single-particle states shown in Fig. 5 are becoming quasi- 
particlc resonances. To describe them properly we have 
solved the HFB equations with scattering type boundary 
conditions. In order to reduce the numerical effort we 
have employed scattering boundary conditions only for 
the quasiparticle states with the energy — A < E < 15 
MeV while for the other states we have used box bound- 
ary conditions. 

With the continuum HFB solutions we can calculate 
how the occupation probability is changing in the region 
around a resonance. This information is provided by the 
quantity 



W(E) = [ 
Jo 



drr* v E (r) , 



(4) 



where v is the lower component of the HFB wave function 
and R is take equal to 20 fm. As an example, in Fig. 6 
we show the values of W(E) corresponding to the states 
3pl/2, 3p3/2 and 2/7/2 in the nucleus 132 Zr. In this 
nucleus the quasiparticle continuum starts at the energy 
E = -A = 0.251MeV. 

Integrating the function W(E) over an energy interval 
in which it has a significant value we can associate to 
each resonance an occupation probability n, i.e.. 



E2 

-- I W{E)dE 



(5) 



The occupation probabilities of the relevant resonant 
states in Zr isotopes are shown in Fig. 7. It can be seen 
that the occupation probabilities corresponding to the 
weakly bound states are increasing progressively when 
passing from 124 Zr to 138 Zr. Thus these states contribute 
significantly to the pairing correlations. 
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FIG. 7: Occupation probabilities for Zr isotopes, calculated 
with SkI4-HFB. 
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FIG. 9: HF (circles) and HFB (squares) neutron radii for Zr 
isotopes. The interaction is SkI4. 




FIG. 8: Contributions of different channels to the to- 

tal neutron density for 124 Zr (top), 132 Zr (middle) and 138 Zr 
(bottom). The interaction is SkI4. 



In order to analyze the structure of the halo, in Fig. 8 
we plot the quantity 



Rijir) 



p(r) 



(6) 



where p and pij are the total density and the density 
corresponding to the channel respectively. From 

Fig. 8 we can clearly see that at large distances the 
dominant contribution to the neutron density is given by 
the p states, which are less confined by the centrifugal 
barrier compared to the other states with higher 
values. This structure of the giant halo obtained by using 
the SkI4-HFB model is very similar to that given by the 
relativistic calculations 0,^2- 



C. The anti-halo effect 

We turn now to the analysis of the so-called " anti-halo 
effect" , which was mainly discussed in light nuclei close to 
the neutron drip line [l5j . This effect is associated with a 
reduction of the neutron radii by the pairing correlations. 
Here, we will show that in fact the pairing correlations 
not only decrease but can also increase the neutron radii 
in nuclei close to the drip line. 

A strong effect of the pairing correlations upon the neu- 
tron radii is expected when close to the Fermi level one 
finds weakly bound s or p states together with states with 
higher angular momenta (weakly bound or resonances). 
As an illustration we present below the case of Zr and Ni 
isotopes. 

In Zr isotopes the pairing correlations affect the neu- 
tron halo through the weakly bound p and / states. In 
Fig. 9 are shown the neutron radii obtained in the HF 
and HFB calculations performed with the force SkI4. We 
observe that in some cases (A = 124, 128, 130, 136) the 
HF radii are larger than the HFB ones while in other 
cases (A = 126, 132, 134) the contrary happens. This 
behaviour can be easily understood if one considers the 
occupancy of the states 3pl/2, 3p3/2 and 2/7/2 in HF 
and HFB calculations shown in Table I and Fig. 7. We 
observe that for A = 124 the state 3p3/2 is half occupied 
while 2/7/2 is empty in HF. When pairing correlations 
are switched on, the occupancy of the state 3p3/2 is re- 
duced and that of the 2/7/2 is enhanced. Since the state 
2/7/2 has a smaller radial extension than the state 3p3/2 
due to the centrifugal barrier, the HFB radius is reduced 
compared to the HF radius. Similar considerations ap- 
ply for A = 128 and 130. For A = 136 the three states 
are completely occupied in HF. The HFB radius is thus 
smaller than the HF one due to the reduction of the oc- 
cupancy of the two p states when pairing is switched on. 

On the other hand, for A = 126 the situation is oppo- 
site: in HF approximation the state 3p3/2 is empty and 
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A 124 126 128 130 132 134 136 

3pl/2 2_ 

3p3/2 2 4 4 2 4 4_ 

2/7/2 4 2 4 8 8 8 



TABLE I: Number of neutrons in the states 3pl/2, 3p3/2 
and 2/7/2 in HF as a function of N in Zr isotopes. 




86 88 90 92 94 96 

A 



FIG. 10: HF (circles) and HFB (squares) neutron radii for Ni 
isotopes. The interaction is SkI4. 



A 86 88 90 92 94 96 98 
3al/2 2 2 2 2 2 2 2 
2ri3/2 2 4 4 4 4 4 
lg7/2\ o|o|o|2|4|6|8 

TABLE II: Number of neutrons in the states 3sl/2, 2d3/2 
and l<?7/2 corresponding to Ni isotopes and HF calculations 
with the SkI4 force. 

which all the states involved in the anti-halo effect are 
bound states. 

Coming back to the neutron radii shown in Fig. 10, we 
can see that, from 86 Ni to 94 Ni the HF radii are larger 
than the HFB ones, the difference being maximum for 
90 Ni. As in the case of Zr isotopes, the effect of the pair- 
ing correlations upon neutron radii ca be simply traced 
back to the occupancy of the HF levels shown in Ta- 
ble II. The most important thing we can notice in Ta- 
ble II is that the state 3sl/2 is fully occupied in HF 
for all the considered isotopes. Thus in all these nuclei 
the pairing force will depopulate the state 3sl/2, which 
has the largest spatial extension compared to the other 
high-/ states shown in Table II. This explains why in Ni 
isotopes the HF radii are systematically larger than the 
HFB radii. 



the state 2/7/2 is partially occupied. Thus, in this case 
pairing enhances the occupation of the state 3j»3/2 and 
reduces that of the state 2/7/2. Consequently, the radius 
becomes larger in HFB than in HF. A similar mechanism 
applies for A = 132 where the state 2/7/2 is completely 
occupied in HF and 3p3/2 only partially occupied. For 
A = 134 these two states become completely occupied 
in HF. The pairing now scatters some neutrons in the 
state 3pl/2 which is empty in HF. The global effect is 
a very weak enhancement of the neutron radius in HFB. 
In conclusion, in Zr isotopes the pairing correlations can 
reduce or increase the neutron radii according to the rel- 
ative occupancy of the weakly bound states close to the 
Fermi level. 

We complete our analysis by considering the Ni iso- 
topes close to the neutron drip line. The anti-halo effect 
in Ni isotopes was discussed already in Ref.0 by using 
the Skyrme force SIII. For this force the neutron drip 
line is extended up to the nucleus 88 Ni. This prediction 
is very different from the RHB calculations, which give 
98 Ni as the last bound nucleus. Interestingly enough, the 
same position of the drip line, i.e., 98 Ni, is predicted by 
the force SkI4. 

The neutron radii of Ni isotopes calculated with the 
SkI4 force are shown in Fig. 10. For the Ni isotopes 
analyzed here the most important single-particle states 
are 3sl/2, 2<i3/2 and lgl/2. The first two states are 
weakly bound while l</7/2 is a resonance state in 86 Ni 
and 88 Ni and becomes weakly bound beyond 88 Ni. This 
level structure is different from the case of Zr isotopes in 



IV. CONCLUSIONS 

In this work we have examined the evolution of the 
nuclear structure of Zr isotopes at large neutron excess, 
with a special attention to the far out region of neutron 
densities. This is motivated by the predictions of the 
RHB approach which indicate the presence of a giant 
neutron halo in these nuclei. We have therefore used a 
different approach, namely the Skyrmc-HFB model. We 
find that, in this isotopic chain the presence or absence 
of giant halos depends essentially on the location of the 
predicted neutron drip line. Thus, there is a strong model 
dependence in this type of study as illustrated by the 
results obtained with the Skyrme forces SLy4 and SkI4. 

For the drip line to be displaced towards heavier iso- 
topes, a necessary condition is that some HF orbitals 
become bound when A increases. In this case additional 
bound neutrons can be accommodated and bound nuclei 
of heavier mass can be formed. An illustration of this 
situation is provided by the Zr isotopes where this nec- 
essary condition is fulfilled by some model like SkI4 but 
not by SLy4. Once the necessary condition is realized, a 
neutron halo may exist if some of the weakly bound HF 
orbitals correspond to low angular momenta (3p3/2 and 
3pl/2 in the case of Zr with SkI4) so that the centrifugal 
barrier is weak enough to let the wave functions extend 
far out. 

Thus, the decisive factor is the HF mean field while 
the pairing correlations play a lesser role. Which part 
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of the mean field governs the position of the drip line is 
not clear, and this point deserves further study. On the 
example of Zr with SkI4 and SLy4 it seems that the main 
role cannot be attributed to the spin-orbit potential but 
rather to some part of the mean field which influences 
the relative positions of lh and 3p orbitals. 

We have also seen that the pairing correlations can 



lead to the anti-halo effect, as it is well known, but it 
can also sometimes produce an enhancement of the halo 
effect. This can be understood by analyzing the occupa- 
tion probabilities of the least bound (lj) orbitals. Finally, 
it can be noted that the anti-halo effect can occur even 
when all active orbitals arc bound. 
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